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Abstract 



Controlling the translational motion of cold atoms using optical lattice potentials is of both 
theoretical and experimental interest. By designing two on-resonance time sequences of 
kicking optical lattice potentials, a novel connection between two paradigms of nonlinear 
mapping systems, i.e., the kicked rotor model and the kicked Harper model, is established. 
In particular, it is shown that Hofstadter's butterfly quasi-energy spectrum in periodically 
driven quantum systems may soon be realized experimentally, with the effective Planck con- 
stant tunable by varying the time delay between two sequences of control fields. Extensions 
of this study are also discussed. The results are intended to open up a new generation of 
cold-atom experiments of quantum nonlinear dynamics. 

Key Words: kicked rotor model, optical lattice, ultracold atoms, Hofstadter's butterfly 
spectrum, kicked Harper model 
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I. INTRODUCTION 

One main objective of the field of quantum control is to use controlled laser-matter 
interaction to explore new aspects of quantum dynamics and enhance our understanding of 
quantum coherence phenomena. Along this direction quantum control ideas and techniques 
are expected to be very useful for quantum simulation studies, i.e., using controlled quantum 
systems to simulate important models. To that end ultracold atoms and molecules offer 
promising opportunities due to experimental advances in Bose-Einstein condensation, the 
great controllability of ultracold systems by laser fields, and the long decoherence time of 
ultracold systems. 

Motivated by our interest in understanding quantum coherence effects and their control 
in classically chaotic systems 3J, our work here focuses on possible cold-atom realizations of 
fundamental models of quantum nonlinear dynamics. One important paradigm of quantum 
chaos is the kicked rotor model (KRM) [4|, |5( whose scaled Hamiltonian can be written as 

2 

H krm = V - + Kcos(q) - nT ), (!) 

n 

where q (g [0,27r)) and p are conjugate coordinate and momentum variables, and T is the 
period of the kicking potential. For a general value of T, the spectrum of H KRM is dis- 
crete and the quantum diffusion of the momentum distribution saturates as a consequence 
of quantum destructive interference despite the unlimited classical diffusion. This has been 
well understood in terms of a one-dimensional Anderson disorder model |6| . Interestingly, if 
the period of the kicking potential T is on resonance with the quantum recurrence time of 
the quantum free rotor dynamics, i.e., T is a rational multiple of the latter, then the asso- 
ciated spectrum will in general consist of continuous bands and ballistic quantum diffusion 
emerges. These features of the quantum KRM have played a significant role in advancing our 
understanding of quantum nonlinear dynamics. They have also motivated the experimental 
realization of the KRM using cold atoms subject to kicking optical lattice potentials Q. 
Indeed, with the joint efforts of about ten laboratories worldwide working on the cold-atom 
realization of the KRM , many important dynamical features of the quantum KRM 

and its variants have been observed. 

Another paradigm of quantum nonlinear dynamics is the kicked Harper model (KHM) 
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Il4 Il5| Using similar notation as the KRM, the KHM Hamiltonian is given by 
Hkhm = (L/T) cos(p) + Kcos(q) ^5(t - nT), (2) 

n 

where L and K are two system parameters. For later use we note that the classical KHM 
map is given by 

p c (n+l) = p c {n) + Ksm[q c (n)], 

q c (n + 1) = q c {n) - L sm[p c {n + 1)], (3) 

where [q c (n) , p c {n)] denote the values of the classical coordinate and momentum right before 
t = nT. The quantum map associated with each period T is given by 

- L cos(p) ■ K cos(g) 

Urhm = e~ l ^^e~ 1 ^^, (4) 

where H is the effective Planck constant in our scaled unit system (hence p = —ihd/dq). For 
the critical case of K = L, the (quasi-energy) spectrum of the quantum map Ukhm is very 
similar to the famous Harper model for studies of two-dimensional electron gases in a strong 



magnetic field [16J. In particular, the spectrum of Ukhm with K = L is a fractal, often called 
the "Hofstadter's butterfly" spectrum lf|. As such, the dynamical properties of the KHM 
are often in sharp contrast to those in the KRM. Indeed, complementing the KRM, the KHM 
has been regarded as another important paradigm of quantum chaos, offering a test bed 
for understanding how a fractal spectrum is manifested in quantum mapping systems and 
how the underlying classical chaos affects the butterfly spectrum. One important question 
thus arises before the research community: How to experimentally realize or experimentally 
simulate the KHM? 



Two previous studies have proposed to use Fermi-surface electrons in pulsed fields [23] 
or a charged particle kicked by a designed field sequence [24] to realize the KHM. However, 
these proposals, unrelated to ongoing cold-atom experiments of quantum chaos, have not 
led to experiments. An intriguing connection between the KHM and the so-called kicked 
harmonic oscillator model [25( might also help realize the KHM, but unfortunately this 
connection is subject to the strong restriction of K = L. By contrast, in a recent work 



17| . we have briefly reported that based on existing cold-atom experiments of the KRM, 
it should be possible to experimentally realize a quantum version of the KHM by making 
use of two sequences of kicking optical lattice potentials (i.e., making use of a double kicked 



rotor system, to be explained later). This unveils for the first time a direct connection 
between the two paradigms of quantum nonlinear dynamics, i.e., the KRM and the KHM. 
In this paper, we present further discussions and important details of our finding. We shall 
emphasize that the effective Planck constant of our realization of the quantum KHM can 
be tuned by varying the time delay between two sequences of control fields, thus offering 
opportunities for understanding the differences and similarities between the quantum and 
classical dynamics. Extensions of our study to a wider class of periodically driven quantum 
systems are also discussed. It is our hope that this contribution can motivate more interest 
in quantum control and quantum simulation of classically chaotic systems. 



II. KICKED HARPER MODEL REALIZED BY ON-RESONANCE DOUBLE- 
KICKED ROTOR MODEL 

Our starting point is a modification of the standard KRM, called a double kicked rotor 
model (DKRM). That is, within each period T, the free evolution of a rotor is interrupted 
twice by external kicking potentials. Such a model has been experimentally realized jl8| and 
has already attracted considerable interests 
above, the DKRM Hamiltonian is given by 



In terms of the scaled variables used 



H=— + K 1 cos(g) 5 ( f ~ nT ) + K 2 cos (<?) Yl S ^ ~ nT -v)- (5) 

n n 

Clearly, within each period T the rotor experiences two kicks at t — nT and t = nT+rj, with 
the two field amplitudes characterized by K\ and K 2 . Hence, there are now two sequences 
of kicking fields with the same period T, and the time delay between the two sequences of 
control fields is given by rj. The associated quantum map Udkrm for a period from nT + 0~ 
to (n + 1)T + 0" is found to be 

Udkrm = e'^-^ cos{q) e'^ cos{q) . (6) 

Note that for a Hilbert space satisfying the periodic boundary condition associated with q — > 
q+2n, which should be the case for a rotor, the momentum eigenvalues can only take integer 
values mutiplied by h. Consider now what happens under the quantum resonance condition, 
i.e., Th = Air. Due to the discreteness of the momentum eigenvalues, one immediately 
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p 2 

obtains e~ lT 2K = 1. Under this resonance condition, Udkrm is reduced to U DKRM , 



u DKRM e c c t 



— i-S 2 - cos(q) — — cos(o) /ryx 

= e 2»e ft e ane « , (7) 

where we have defined the rescaled momentum p = rjp and the rescaled kicking amplitudes 
K\ = r]Ki, K 2 = f]K 2 . In terms of this rescaled momentum operator p, the effective Planck 
constant evidently becomes % = r]%. Certainly, in a cold-atom realization of this system, one 
should not forget that the cold atoms are actually moving in a flat space rather than in a 
compact angular space. Hence, the quantum resonance condition is relevant only when the 
initial quantum state is prepared in a state closely resembling a momentum eigenstate. This 
is already well within reach of today's experiments. For example, two recent experiments 
0, [lo| realized a KRM on quantum resonance, using a delocalized Bose-Einstein condensate 
that generates appropriate initial states. 

Using the rescaled variables defined above and the associated effective Planck constant 
h, the quantum map in Eq. ([7]) can now be interpreted in a straightforward manner. Specif- 
ically, within each period T, the system is first subject to one kick, followed by a free 
evolution interval; then the system is kicked a second time, followed by a second interval 
of free evolution, with the free Hamiltonian for the second free evolution interval given by 
H{ T ee = — P 2 /2. This interpretation makes it clear that the on-resonance double kicked rotor 



model realizes a modified kicked rotor we proposed in Ref . [2l| . Equation ([7]) also indicates 
that the time delay rj between the two sequences of the kicking fields offers a convenient 
means to vary h = rjh. Hence, the effective Planck constant h of the quantum map obtained 
above can be easily tuned, so long as we keep K\ and K 2 constant. This offers a promising 
opportunity for studies of the quantum-classical correspondence associated with the map 

jjr 

u DKRM- 

Consider now a well-defined classical limit of the map U r RDKR , i.e. the 77 — >• limit or 
equivalently the h — > limit. Note that this special classical limit is based on the above on- 
resonance DKRM, and is hence unrelated to the conventional classical version of a DKRM 
(the conventional classical version of a DKRM is reached by letting % — > with fixed Ki, 
K 2 , T, and 77). Based on the above simple interpretation of U DKRM , let us now consider the 
classical analogs of the quantum observables p and q right before t = nT. These classical 
quantities will be denoted p c (n) and q c (n), representing the values of classical canonical 



variables right before t = nT. Then, after the combined action of the first kick at t — nT 
and the first free evolution interval, p c {n) and q c {n) evolve to p' c {n) and q' c (n), with 

p' c (n) = p c (n) +^isin[g c (n)], 

q' c ( n ) = Qc(n)+p' c (n). (8) 

Similarly, after the combined action of the second kick and the second free evolution interval, 
the total elapsed time is T and p' c {n) and q' c (n) evolve to p c {n + 1) and q c {n + 1), 

p c (n + 1) = p' c (n) + K 2 sm(q' c ), 

q c (n + l) = q' c (n)-p c (n + l). (9) 

The overall classical map in the limit of % — > 0, for the period between nT + 0~ and 
{n + 1)T + 0~, can be found by considering the above two steps together, i.e., 

p c (n + 1) = p c (n) + K 2 sin q c (n) + p c (n) + K x sin[g c (n)] + K x sin[g c (n)], 

q c (n + 1) = q c (n) - K 2 sin q c (n) + p c (n) + K x sm[q c (n)} . (10) 

At first glance this map seems rather complicated, but if we make a classical canonical 
transformation 

(Qc,Pc) -> (Qc = q c ,P c = Pc + q c ), (11) 
then the map of Eq. fllUI) assumes a much simpler form, 

P c {n + 1) = P c {n) + Kt sm[Q c {n)\, 

Q c {n + 1) = Q c {n) - K 2 sin[P c (n + 1)]. (12) 

Remarkably, with the substitutions K\ — > K and K 2 — > L, the map we obtained in Eq. 
f[T2"j) is identical with the classical KHM map in Eq. ([3])! That is, for a DKRM under the 
resonance condition Th = 47r, the quantum map U 1 DKRM is simply a quantum version of the 
KHM, insofar as its % — > classical limit is equivalent to a classical kicked Harper model. 
This surprising finding exposes a direct connection between the KHM and a kicked rotor 
system for the first time. Because the DKRM was experimentally realized a few years ago, 
the finding here indicates that our quantum version of the KHM can be realized by slightly 
modifying the previous experiments [3]. Considering the many fascinating dynamical fea- 
tures of the KHM, we hope that our work will be able to motivate many new cold-atom 
experiments of quantum nonlinear dynamics. 



o!o 0.5 j!o 1.5 Z0 



FIG. 1: Quasi-energy spectrum of the cold- atom realization of a quantum version of the kicked 
Harper model, achieved by considering a double kicked rotor model under the main quantum 
resonance condition TTi = 4-7T. For panels (a), (b), and (c), K\/Ti = K2/T1 = 2.0, 3.0, and 4.0, 
respectively. The Hofstadter's butterfly spectrum seen here is almost indistinguishable from that 
calculated from the standard kicked Harper model. 

III. DETAILED RESULTS 

A. Hofstadter's butterfly spectrum 

In Fig. 1 we show the quasi-energy spectrum of U r DKRM as a function the effective Planck 
constant h, for three values of Ki/h = K 2 /h. The famous Hofstadter's butterfly structure 
can be clearly seen. Indeed, we have compared the spectrum of U r DKRM with that of the 
Ukhm hi Eq. (jlj), for K\ = K = K2 = L and for the same value of the effective Planck 
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constant. The result is that no difference in the Hofstadter's butterflies can be seen by 
the naked eye. For other parameters (K\ = K ^ K 2 = L), the remarkable similarity is 
also observed (not shown). The fractal property of the spectrum of U r DKRM has also been 



checked carefully. For example 
generalized fractal dimension 



or K x = K 2 = 1, h = 2tt/(1 + a), and a = (y/E + l)/2, the 
171 ] of the spectrum is found to be Dq w 0.5. This is identical 
to that of the KHM, with K = L = 1 and the same value of the effective Planck constant. 

The results in Fig. 1 bring up an interesting question. Is the butterfly spectrum of U r DKRM 
mathematically the same as that of the KHM? That is, does there always exist a unitary 
transformation to transform Ukhm in Eq. PJ) to U r DKRM l The answer is no. Even though a 
canonical transformation of the classical limit of U r DKRM is equivalent to the classical limit of 
Ukhm, their quantum spectra can still be different due to the periodic boundary condition 
associated with their Hilbert spaces. Qualitatively, this is because a classical canonical 
transformation does not necessarily have a unitary transformation analog in the quantum 
case. In other words, a particular quantization strategy can induce differences between the 
butterfly spectra of U r DKRM and that of Ukhm- 

Especially, we have shown that for fixed K/% and L/h, the spectrum of Ukhm has a 
period 2n in h, and is reflection symmetric about h — it. By contrast, the spectrum of 
U r DKRM (for fixed Ki/H and K 2 /h) has a period An in h, and is reflection symmetric about 
h = 2ir instead. Also significant, it can be proved that the spectrum of U r DKRM is invariant 
if we swap K\ and K 2 . 

As an example of spectral differences, we find that so long as K\ 7^ K 2 , the spectrum 
of U r DKRM is in general continuous. This is markedly different from the KHM, where the 
spectrum can change from being continuous to being discrete if we swap the values of K 
and L [13I, Q- Other subtle spectral differences between U r DKRM and Ukhm, especially 
those regarding the butterfly spectrum's sub-band widths in cases of ft = h — 2nr/s (r and 
s being integers), have also been studied [^j but will not be discussed here. 



B. Quantum Diffusion Dynamics 

Our cold-atom proposal for realizing one quantum version of the KHM has presented ex- 
citing opportunities for experimental studies of quantum diffusion dynamics with a butterfly 
spectrum. Though the fractal butterfly spectrum seen above is not a direct experimental 
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FIG. 2: Time dependence of the momentum variance crj(i) in the cold-atom realization of the 
kicked Harper model. The main features range from localization, subdiffusion to almost ballistic 
diffusion. In the dotted line case h = 67r/19; for other three cases h = 1. The initial state is taken 
as a zero momentum state, which can be approximately realized by use of a dilute Bose-Einstein 
condensate with very large coherence length. 

observable, the fractal properties of the spectrum do govern the quantum dynamics 0. A 
number of previous studies have contributed to understanding the manifestations of a fractal 
spectrum in the quantum dynamics. 

To motivate experimental work based on our finding we present here some remarkable 
numerical results of the quantum diffusion dynamics associated with U r DKRM . We choose 
to examine the time dependence of the momentum variance of a time-evolving wavepacket 
\(fi(t)). This variance is given by erj(t) = h ^2i\(^(t)\l)\ 2 (l — lo) 2 -, where {|/)} are the 
eigenstates of p with the eignevalue lh, \lo) is the initial state, and Iq is the time-independent 
mean value of p/h. Note that crj(£) is a quantity easily measurable in current cold-atom 
experiments of the KRM. In our numerical calculations we assume Z = but we have 
checked that the results presented below do not depend on the initial condition we choose. 

First of all, for both K\ <C K 2 and K\ ^> K 2 , the dynamics of U r DKRM typically displays 
almost ballistic diffusion with crj(t) ~ t 2 . One such example is represented by the dot-dashed 
line in Fig. 2. It is also computationally found that if we swap K\ and K 2 , then the time 
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dependence of erj(t) does not change. This feature is consistent with our previous discussion 
and does not exist in the dynamics of Ukhm- 

Besides, for K\ = K 2 and for a generic value of h, the quantum diffusion dynamics is 
expected to be drastically different because the spectrum is a fractal (see Fig. 1). In Fig. 
2 we show one such computational example with K\ = K 2 = 3.9 and h — 1 (solid line). It 
is seen that the quantum diffusion dynamics in this case is much slower than the previous 
ballistic diffusion case. To further understand this case we follow the dynamics over a time 
scale of t — 10 4 T and then fit cx|(t) by a power law. This fitting yields cx?(t) ~ t a with 
the quantum diffusion exponent a ~ 0.82, indicating an anomalous diffusion. Applying an 
existing theory [22| that relates the quantum diffusion exponent a to the fractal dimension 
of the spectrum, we infer that in the case of a ~ 0.82, the Hausdorff-dimension of the fractal 
spectrum is Dh ~ 0.41. Interestingly, if we slightly mismatch K\ and K 2 (dashed line in 
Fig. 2), then it is found that after a transient time of about 10 2 T, cr|(t) will display ballistic 
diffusion like in the first example. As such, by merely slightly tuning the relative strengths 
of the two sequences of the control fields without even changing their average strength 
(K\ + K 2 )/2, one may generate many remarkable and qualitatively different features in the 
quantum diffusion dynamics. 

Let us finally discuss the result of the dotted line in Fig. 2. In this case, h = 67r/19, 
K\ = K 2 = 1.8, and aj(t) is seen to saturate at very small values. Computationally, this 
saturation behavior is found to persist for times larger than 10 6 T. Hence the quantum dif- 
fusion dynamics in this case constitutes an example of strong localization. This localization 
behavior is related to the fact that the spectrum sub-band width of the butterfly of U r DKRM 



may become vanishingly small as K\ = K 2 decreases 27|. We believe that the localiza- 



tion observed here represents a novel type of dynamical localization in quantum nonlinear 



dynamics. Our ongoing theoretical work will soon offer more insights into this [27| . Experi- 
mentally speaking, observing such kind of localization behavior would be a strong indication 
that our cold-atom version of the KHM has been cleanly realized. 

C. Extension to Double Kicked Rotor Systems on Higher-Order Resonances 

So far, we have shown that by considering a double kicked rotor under the resonance con- 
dition Th = 47r, a cold-atom realization of the KHM is within reach of today's experiments. 
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FIG. 3: Quasi-energy spectrum of a double kicked rotor system under the quantum anti-resonance 
condition TTi = 2ir. In panels (a), (b) and (c), K\j% = K^/fi = 1.2, vr/2, and 2.0, respectively. 
In (b), three main branches of the spectrum pattern begin to touch each other. The complex 
spectrum pattern as a function of U is called "generalized Hofstadter's butterfly" in the text. 

It now becomes interesting to ask whether Hofstadter's butterfly spectrum also exists for a 
double kicked rotor system under higher order quantum resonances, i.e., under the condition 
Th = Airu/fi, with v and [i being integers and with \x > 1. Our preliminary study suggests 
that this is the case, thus opening up studies of quantum diffusion dynamics in a wider class 
of quantum mapping systems. 

As an example, we present in Fig. 3 the spectrum of the DKRM at the so-called quantum 
anti- resonance, with v = 1, /i = 2. For Ki/H = K 2 /h = 0, the quantum map operator 
Udkrm in Eq. (jEJ) reduces to U = e ~ mp2 / %2 , yielding only two nonequivalent values of 
the quasi-energy, i.e., and ±n. As Ki/h = K 2 /h increases, the complex structure of 
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the spectrum "grows" around the points and ±7r. Due to this growth, the spectrum for 
the case of Ki/K = K 2 /h = 1.2 is already quite complicated, as seen in Fig. 3(a). When 
Ki/h = K 2 /h reaches n/2 [see Fig. 3(b)], three main branches of the spectrum begin to touch 
one another, yielding an overall pattern that can be regarded as a generalized Hofstadter's 
butterfly spectrum. With the values of K\/% = K2/T1 increased even further, the generalized 
Hofstadter's "butterfly" pattern gets larger and possesses more fine structures. Unaware of 
any study with similar results, we think that the generalized "butterfly" spectrum found here 
and those associated with other higher-order quantum resonances will stimulate considerable 
theoretical work. The results further strengthen the view that using controlled laser-matter 
interactions, novel quantum dynamics models may be generated and explored. 



IV. CONCLUSION 



To conclude, we have shown that by considering a cold-atom realization of an on- 
resonance double kicked rotor model we can realize a quantum version of the kicked Harper 
model from a kicked rotor system. To the naked eye, the butterfly spectrum we obtain 
from our version of a kicked Harper model is almost indistinguishable from that of the stan- 
dard kicked Harper model. We have also stressed that the effective Planck constant of the 
quantum kicked Harper model realized here can be easily tuned by varying the time delay 
between two sequences of control fields. Extending our considerations to double kicked rotor 
systems under higher-order quantum resonances, we have shown that an entirely new class 
of quantum mapping systems with generalized Hofstadter's butterfly spectra can be studied, 
both theoretically and experimentally. These results present many new opportunities in the 
studies of quantum nonlinear dynamics in external control fields. Indeed, with controlled 
interactions between laser fields and cold atoms, there are still so much to see and explore. 
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